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We present an investigation of various gauge invariant definitions of the qq Bethe- 
Salpeter (BS) amplitude for mesons in lattice QCD, and compare them to the Coulomb 
and Landau gauge BS amplitudes. We show that the gauge invariant BS amplitude is 
considerably broadened by the use of "fat" gauge links (constructed by smearing the links 
of the original lattice). A qualitative demonstration of the Lorentz contraction of the 
pion wavefunction at non-zero momentum is given. We also calculate density-density 
correlations and discuss the limitations in extracting the charge radius of the pion from 
these observables. Lastly, the polarization dependence of the BS amplitude for the p meson 
is exhibited, and we extract the relative admixture of / = and 1 = 2 states showing that 
simple hydrogen-like wavefunctions provide a good fit. 
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1. Introduction — Bethe-Salpeter amplitudes 

In the study of bound states in relativistic quantum field theory, commonly used 
tools are the Bethe-Salpeter (BS) amplitudes. These give the probability amplitude for 
finding some specified arrangement of fundamental components within the bound state. 
In this paper we study three different types of equal-time BS amplitudes for the pi and rho 
mesons and discuss whether one can extract an experimentally measured quantity such as 
the charge radius from such probes. 

A gauge invariant definition of BSA for, say, a pion of momentum p is 

A^x) = (0K0)75A^(0,x)u(x)|7r(p)) (1.1) 

where A^(0,x) is a path-ordered product of gauge links that joins points x and and 
makes the amplitude gauge invariant. This amplitude is given by the following ratio of 
2-point correlators: 

^ (0|40; t)yM{Q, X- t)u{x- 1) u{y- Q)^,d{y- 0)|0) 

" (0|(/(0;t)75u(0;t) u(y;0)75(/(y;0)|0) 

M. may be defined in a number of ways. In this paper we use the APE smearing method [1] 
to construct "fat" gauge links and thus obtain a nonlocal average of paths for M.. We also 
consider the Coulomb gauge {C-^) and Landau gauge {C-^) BS amplitudes. These are 
obtained by transforming the quark propagators to Coulomb (Landau) gauge, and then 
calculating the ratio given in Eq. (1.2) with M. = 1, i.e. without including the links. The 
pion decay constant /^r is related to ^7r(0) and has already been calculated for this data 
set in Ref. [2]. Here we are only interested in the behavior of A-k{x) as a function of x so 
it is expedient to use the normalization ^7r(0,t) = 1. 

Calculations of BS amplitudes in lattice QCD began with the 1985 study of Velikson 
and Weingarten [3]. Other calculations include the work of Chu eA al. [4], Hecht and 
DeGrand [5], Duncan eA al. [6], and Kilcup [7]. Most of this work has used the Coulomb 
or Landau gauge, though there is some work on simple gauge invariant prescriptions [4] [5]. 
A preliminary version of the results presented in this paper have been reported earlier in 
Ref. [8]. Another similar probe of the hadronic structure are density-density correlations, 
as suggested in Refs. [9] [10] [4]. We discuss this method in section 5. 

For gauge theories, the definition of BS amplitudes is complicated by the problem 
of gauge dependence. In QED the gauge dependence can be calculated and is typically 
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0{aem) [11]- Therefore one is free to choose a gauge that makes calculations simple. This 
is why it is common to choose Coulomb gauge for QED calculations as it has the useful 
property of regulating infra-red divergences [12]. One does not expect such a simplification 
in a strongly coupled theory like QCD, instead to measure the distribution of constituents 
inside a hadron one needs a probe that takes into account the chromo-electric flux and qq 
pairs in addition to the valence quarks and antiquarks. 

In the case of mesons we imagine that the valence q and q are effectively conflned by 
a flux tube. (In this picture the effect of additional qq pairs from the vacuum is to simply 
decrease the value of the string tension.) Since it is known from experiments that the 
glue carries a signiflcant fraction of the hadron's momentum, different deflnitions of BSA 
will give very different answers depending on the overlap of the probe with the physical 
state. For this reason we characterize the "goodness" of a probe by the magnitude of the 
amplitude. 

The virtue of Coulomb and Landau gauges is that they make color flelds smooth by 
reducing the gauge fluctuations, so the flxed gauge BSA measures that part of the meson 
wavefunction corresponding to a smooth distribution of chromo-electric flux between the 
qq. The principal disadvantage of calculations in these gauges is the poorly understood 
effect of Gribov copies. A review of the Gribov ambiguity in lattice calculations can be 
found in Ref. [13]. On the other hand it is obvious from Eq. (1.1) that in gauge invariant 
deflnitions each different choice of the connection Ai will give a different BSA. By using 
smeared links to construct the path joining qq we are using a certain linear combination of 
paths on the original lattice. There is, however, little control over the relative weights of 
the different paths, and in this sense the smearing process is not optimal. The usefulness of 
the smearing process lies in the fact that the overall thickness of the flux tube in the probe 
is controlled by the number of smearing steps and that the method is computationally 
simple. 

It should be emphasized that while the BS amplitudes are simple to measure, there 
is, for light quarks, no demonstrated connection between a hadron's size extracted from 
them and the experimentally measured charge radius. What BS amplitudes do provide 
is a qualitative understanding of hadronic structure. This knowledge is useful for lattice 
studies in designing better probes, i.e. those that have a larger projection onto a given 
state. This technical improvement is crucial for improving the calculation of the spectrum 
and of matrix elements between hadronic states. 
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This paper is organized as follows. In Section 3 we present results for the gauge 
invariant BS amplitude as a function of the smearing size of the links and show that 
the amplitude becomes broader on using "fat" links to connect the qq. We also give 
a qualitative demonstration of the Lorentz contraction of the wavefunction at non-zero 
momentum. A comparison between Coulomb gauge, Landau gauge and gauge-invariant 
amplitudes is made in Section 4. Calculation of density-density correlations is presented 
in Section 5 along with a discussion of why the charge radius measured in experiments 
is different from similar quantities that can be extracted from lattice data. In Section 6 
we show that the BS amplitude for the p meson depends on the polarization axis. We 
extract the relative admixture of / = and 1 = 2 states from this data and find that simple 
hydrogen-like wavefunctions provide a good fit. Final conclusions are given in section 7. 

2. Details of the lattices and propagators 

The calculation is done using 35 background gauge configurations of size 16^ X 40 
at fi = 6.0. The Wilson action quark propagators are calculated on doubled lattices 
(16^ X 40 — > 16^ X 80) using Wuppertal sources. The quark masses used are k = 0.154 
and 0.155, corresponding to pions of mass 700 and 560 MeV respectively. Further details 
of lattice generation and propagator inversion are given in Ref. [2]. The main limitation 
of this calculation is that the spatial lattice size, L = 16, is not large enough to reliably 
extract the asymptotic behavior for some BS amplitudes and density-density correlations. 
We shall address this issue at appropriate places in the analysis. 

We use a stochastic over-relaxed {u = 2) algorithm for gauge fixing to either Coulomb 
or Landau gauge. In the iterative process, the new gauge transformation at a given site is 
set to either Gnew = AG * GqU or Gnew = AG * AG * GoU with equal probability. The 
change AG is the SU(3) matrix that maximizes Tr (AGS^) where E is the sum of links 
originating from that site. The convergence criteria used in both cases is (3 — Tr AG)/ 3 < 
10~^. Operationally, we first fix the lattice to Landau gauge and then fix to Coulomb 
gauge as this saves computer time. 
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3. Gauge invariant BS amplitude with smearing 

In an earlier calculation Chu et al. [4] investigated the simplest version of the gauge- 
invariant BS amplitude; for A^(0,x), they used the straight line path between points that 
lie along one of the lattice axes. They found that {r'^)n measured from the simplest gauge 
invariant BS amplitude is smaller than that obtained in either Coulomb or Landau gauge 
and that the latter estimates are 0.3-0.5 of the charge radius measured in experiments. 
Similar results have also been obtained by Hecht and DeGrand [5]. 

The poor quality of these earlier gauge invariant results is due to having used a 
straight line path to join the quark and anti-quark: such a probe has a poor overlap with 
the hadronic state because there is only a small probability amplitude for the gluon field 
to be so localized. The overlap can be improved significantly by using a smeared gluon 
field, and the resulting gauge-invariant BS amplitude is markedly enhanced. 

To construct "fat" links we use the APE smearing method, which was first introduced 
to enhance the signal in glueball calculations [1]. In this method each link in the spatial 
direction i is replaced by the sum 

u\'^\x,x + i) = v(^U,{x,x + i) + ^S,{x,x + i)\ (3.1) 
^ i=i ^ 

where Si are the four spatial staples shared by the link ?7i, and the symbol V implies 
that the sum is projected back on to the group SU(3). As a result of this projection 
the normalization of the BS amplitude on different smeared levels is the same, so we can 
directly compare the results. The smeared link represents locally averaged gauge fields and 
has the same gauge freedom as the original link, i.e. under Ui(x) — )• g{x)Ui{x)g{x -\- i) the 
smeared link also transforms as U^^\x) — )• g{x)U^^\x)g{x -\- i). Note that the new lattice 
is of the same size as the original. One can iterate this smearing step as many times as 
necessary, using the effective fields at step (n — 1) to produce still "fatter" fields at step 
n; for example in the second smearing step the right hand side of Eq. (3.1) is constructed 
from smeared links produced in step one. A straight line path between the q and the g 
that is made up of smeared links is in fact an average over a large number of paths on 
the original lattice. We specify the smearing level by a superscript on A, which will be 
— 6 corresponding to the original links and six levels of smearing. We do not consider it 
appropriate to go beyond 6 levels of smearing on a lattice of size 16 with periodic boundary 
conditions. 
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It is worth mentioning that we had hoped to show that the chromoelectric flux is 
more spread out at smaU separations and gets conflned to a flux tube at large separation 
by contrasting data at the diff'erent smearing levels. Unfortunately, the present data do 
not show any such eff'ect. 

We flrst study the behavior of the BS amplitude as a function of the time-slice t 
from the source in order to determine how fast it converges. In Fig. la we show data 
for Att for t = 2,5, 10 and 15, and similar data for the rho in Fig. lb. The qualitative 
shape of the BS amplitude is similar for all t, though quantitatively it gets signiflcantly 
broader, reaching an asymptotic value by about t = 15 for both the rho and the pion. (The 
difference between t = 10 and t = 15 is signiflcant for the pion but not for the rho.) This 
separation is somewhat larger than t = 10 at which we flnd the onset of the plateau in the 
effective mass plots (which is taken as evidence that the correlator is dominated by the 
lightest state) as shown in Ref. [2] using the same set of lattices. The errors in the data are 
independent of t for the pion and increase with t for the rho. In Figs. la,b, we show the 
jackknife errors at all x at t = 2 and for a;>10att = 15as typical examples. We consider 
the t = 15 data best with respect to extracting the ground state and the statistical signal, 
and quote all subsequent results for this value of t. 

In Fig. 2a, b we show data for A^^\x, t = 15) for the pion and the rho as a function of 
the smearing level i. We flnd that the BS amplitude falls off less rapidly as i is increased. 
The statistical errors are similar for smearing levels 1 — 6 and for clarity we only show these 
for i = 1. The data show a rough convergence by z = 6, however, this is speciflc to the 
particular smearing method we have used and even that needs to be conflrmed on a larger 
lattice and with further smearing levels. Given the success of this simple minded "fat" Ai 
in increasing the overlap with the pion state, we feel it is worthwhile to investigate other 
more physically motivated probes in the future. For present we take the results at z = 6 
as our best estimate for a gauge-invariant BS amplitude. 

One advantage of measuring the gauge-invariant BS amplitude is that we can follow 
the signal out to x = 15 on a L = 16 lattice. The data for the pion (Fig. 3) show 
that for Att the large x behavior is flt well by an exponential for 6 < x < 12, while 
such a behavior is hard to extract from A^-^ K The data for x > 12 show some curvature 
which could be due to neglecting power law corrections (discussed later) or due to flnite 
size effects induced by smearing. Our best parameterization of the asymptotic behavior is 
A';^'> ~ ^-o.30ii)x ^(6) _ g-o.29(i)x at /c = 0.154 and 0.155 respectively. We contrast 
this rate of fall-off with that obtained from density-density correlations and the pion mass 
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in section 5. Data for the rho are analyzed in Section 6 after we formulate the dependence 
of the BS amplitude on the polarization axis. 

A measure of lattice discretization effects can be obtained by calculating these ampli- 
tudes at non-zero momentum with the separation x taken to be parallel or perpendicular 
to the direction of p. If the BS amplitude of a meson at rest is characterized by y, z), 
then under a boost in, say, the z direction the measured amplitude should be given by 
tjj^x^y^jz) where 7 = [1.47, 1.66] is the Lorentz contraction factor at the two values of k. 
We show the data for the pion at k, = 0.155 in Fig. 4. Jackknife errors are shown only for 
the cases x\\p and p = for clarity. The data for the case x_Lp are consistent with those 
for p = as expected. The fall-off for the case x\\p is much faster and the signal extends 
only to x = 6 at either value of k. This range is not long enough to extract an effective 
7; all we can say at present is that qualitatively it is consistent with the expected Lorentz 
contraction. 

4. Comparison between gauge invariant, Landau, and Coulomb gauge BS 
amplitudes 

We have measured the Coulomb and Landau gauge amplitudes for the following rel- 
ative separations: the anti-quark's position is varied in a cube of size L/2 with respect 
to the position of the quark which is taken to be at (0,0,0). For each of these relative 
separations we sum the quark's position over the time slice to produce a zero-momentum 
state. We find that the data for x > 6 along any of the axes show that there is a signifi- 
cant contamination from wrap-around effects due to periodic boundary conditions. These 
effects can be taken into account by subtracting the contributions of all the mirror points. 
In this paper we only use the Landau and Coulomb gauge data for comparison; for this 
purpose it suffices to choose data with \xi < 6| to avoid significant finite size effects. 

The data show that A^^^ is approximately the same as C and slightly broader than 
C. On the other hand we find £ > C > consistent with the earlier results of Ref. [5]. 
These two features are illustrated in Fig. 5a by data for Cjr, A^-^ ^ and A^-^ K Similarly, in 
Fig. 5b we show Ct^-, A^^^ and Ai^\ The slight increase in Ct^- and at large x is most 
likely due to mirror contributions that arise due to the use of periodic boundary conditions. 

Using six levels of smearing on a 16^ lattice with periodic boundary conditions pro- 
duces a linear combination of paths in Ai that corresponds to a fairly smooth distribution 
of gauge fields across the complete time-slice. It is therefore not surprising that C, C and 
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^(6) gjyg very similar results. It is still advantageous to work with gauge invariant BS 
amplitudes as one can probe larger separations (up to x = 15) without having to worry 
about contributions of mirror sources, and even more importantly because there are no 
ambiguities associated with Gribov copies. 



5. Comparison with density-density correlations 

The density-density correlation for a charged pion is defined by 

p"^(x)= I d'y{Q\J^{T)u{y,t)VMM'd{x + M^pd{x (5.1) 



where J^r is a pion source and u(y^t)Tau(y^t) probes the u quark at space-time point 
iy^t). In particular, the correlation between charge density operators, Fq, = = 74, is a 
measure of the charge distribution as a function of separation between quarks [9] [10]. The 
charge radius measured in experiments is defined as 



dF{q^) 



(5.2) 

g2 = o 



where F{q^) is the pion form factor. It was shown in Ref. [14] that the Fourier transform 
of p^^{x) is 

p'\q) = J2 {<P = = OMq)) {n(q)\.h(x = 0)|7r(p = 0)) (5.3) 



n 



where n(q) is a complete set of states. Assuming that only n = tt intermediate state 
contributes, one can show that 



dp'\q) 



dq^ 



dq^ 



(5.4) 

g2 = o 



g2 = 

However, due to the contributions of, say, excited P-wave states there is no simple relation 
between the charge radius and p'^'^ [15] [16]. It is not known whether these corrections 
enhance or decrease the rms radius calculated from p'^'^ with respect to the experimental 
value of the charge radius. Using the data for p'^'^ presented below we find a value that is 
too large, but we remind the reader that the lattice calculation has a number of systematic 
errors that are not accounted for. 

We calculate p'^'^(x) using Wuppertal source quark propagators. The meson sources 
are at time T = and T = 40 corresponding to forward and backward propagation on 
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our lattice. Even though we have measured p'^'^(x) for both the pion and the rho, we only 
present results for the pion as the signal for the rho is too noisy. To improve the signal 
we present results for the ratio TZ"^{x) = p"^{x)/ p"^{0). By dividing by the correlation 
function at zero separation we have also removed uncertainties due to renormalization of 
operators on the lattice and can concentrate on effects which depend on the separation 
between the two density operators. 

The Wick contractions for the density-density correlation yield the three types of 
diagrams shown in Fig. 6. Fig. 6a shows both insertions on valence quarks, Fig. 6b shows 
the correlation between insertion on a sea quark loop and a valence quark within the 
pion, and Fig. 6c shows the correlation between two sea quarks within the pion. All three 
diagrams can be estimated in the quenched approximation, however, in order to include 
diagrams in Figs. 6b, c we require additional calculation of quark propagators from every 
spatial point. This is beyond the scope of the present calculation. Furthermore, each 
of these three diagrams is modified by sea quark contributions; these remain a source of 
systematic error that we cannot yet estimate. 

Hueristically, we expect the effects of diagrams shown in Figs. 6b, c to be suppressed. 
The loop tpjitp can interact through gluon exchange with either a valence quark in the 
pion or with a second loop. At short distances, and when the pion propagates as shown 
schematically in Fig. 7a (with the density insertions replaced by disconnected loops), these 
interactions are suppressed by powers of a^. On the other hand for cases corresponding to 
propagation shown in Fig. 6c the exchange is suppressed by e~'"s^/e~'"'"^ where nig is the 
mass (~ 1 GeV) of the mixed glueball meson state with the appropriate quantum numbers. 
Needless to say, it is primarily for computational reasons that we restrict ourselves to the 
connected diagram of Fig. 6a which can be labeled as the valence density correlation in 
the quenched approximation. 

The asymptotic behavior of the connected density correlation is governed by ampli- 
tudes which fall most slowly as a function of spatial separation x. In our setup, the time 
separation between the source (sink) and the density probes is > 10a, so, based on our 
analysis of the spectrum and the BSA, we are probing the lowest state in the pion channel. 
This t separation is also larger than the maximum spatial separation allowed between the 
density probes, i.e. Xi = 8, so we neglect transient effects due to the spatial position and 
character of the pion source (sink). The leading contributions are then described schemati- 
cally by the two diagrams shown in Figs. 7a, b. The diagram in Fig. 7a shows that the pion 
can propagate by preserving its identity all the way, so 'R.'^'^{x) should fall off as x~"' e~"^^^ 
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as x — )• oo. The factor comes from the superposition of momentum excitations in the 
propagation of any state created by a locahzed source, and in three dimensions the asymp- 
totic form is x~^^'^. In Fig. 7b only one of the quarks propagates in the spatial direction, 
effectively forming a flavor singlet meson with a 7^ point source (under a rotation of axis 
74 — ^ tO- Thus the fall off is expected to go as x~^^'^e~"^f^ for large x. So for pions the 
scenario in Fig. 7a produces the dominant behavior at large distances, while for the rho 
both diagrams fall off at the same exponential rate, i.e. e~"^f^ . 

We measure p'^^(x) at all x in the flrst Brillouin zone. Nevertheless, we flnd that it 
is necessary to account for the leakage between adjacent Brillouin zones due to the spatial 
periodicity of the lattice [17] [18]. The lattice correlation, for a cubic lattice of side L, 
denoted by p(x)l^ is an inflnite sum of all the mirror images 
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p{^)l = 2^ p{x + Ln)oo ■ (5.5) 

ni ,112 ,113 = 00 

In practice we truncate the sum and consider contributions only from \xi — ni \ < L. This 
truncation is justifled since the density correlation falls exponentially at large separations 
with a decay rate equal to the pion mass, thus the contribution from the nearest images is 
roughly e~"^^^ e~^'^ with a maximum degeneracy factor of 8. To perform the correction 
for the periodic images we parametrize p'^^(x) ~ e~"^^^ (we ignore the additional power 
law factor as our data are not sensitive to it) self-consistently until it gives the best flt to 
the data. Further details of this method for image corrections are given in a forthcoming 
paper[17]. 

The raw data for the density correlation and the Bethe-Salpeter amplitude are shown 
in Fig. 8. To interpret these results we flrst deflne and extract an rms radius from each 
quantity and then discuss the physical meanings of these. Assuming Eq. (5.4), the charge 
radius is given in terms of the density-density correlation by 

(rL) = I f d\x'p'\x) (5.6) 

where x is the relative separation of the d and u. Similarly, we deflne an rms radius from 
the BSA amplitude as 

, . . _ l !d^xx\A^){x)f 

^'^^^^ " w Jd^x(A(^)(x)f ^^-^^ 
where the factor w translates the relative separation to separation from the center of mass. 
The two limiting cases are w = A for a two body system with degenerate q and g, and 
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w = 2 if the q and q move independently about the center of mass. The real situation is 
somewhere in between, and at present we do not have a way of estimating it [15]. We shall 
calculate vbsa assuming w = 4 and comment on this uncertainty later. 

To calculate these rms radii from our data requires that we make an ansatz for 
the asymptotic fall-off. For the density-density correlations we use both e~"^^^ and 
x~^^'^e~"^^^ where m^r is the pion mass extracted from the 2-point function. We regard the 
difference between the results from the two ansatze a measure of the uncertainty induced 
by this process. For the BSA we simply use the fits shown in Fig. 3. 

The final results are given in Table 1, and these are to be compared with the ex- 
perimental value = 0.636 it 0.037 fm [19]. We also include for comparison results on 
16"^ lattices at /3 = 5.7. A description of this data set is given in Ref. [18], and we have 
reanalyzed it as described in this paper. A comparison of rdd obtained at the two values 
of /3 is particularly significant because the pion mass is roughly the same in physical units. 
This is true at each of the two values of the quark mass. The data show a significant 
variation with the quark mass and extrapolating to the physical pion mass will give an 
rdd that is significantly larger than the experimental value at both values of /3. Chiral 
perturbation theory predicts that r^r has a logarithmic singularity in the chiral limit [20]. 
In the quenched approximation, however, the pion loop that gives rise to this singularity 
is absent, so we expect a finite value for the rms radius in the chiral limit. The data 
also show that the deviation is significantly reduced on going from /3 = 5.7 — > 6.0; to get 
the continuum result one needs to work at still weaker coupling or to use an improved 
action. Also, the effects of the disconnected graphs, i.e. Figs. 6b, c, and quenching need to 
be investigated before we can quantify the effects of the neglected (n ^ vr) states in Eqs. 
(5.3) (5.4). 

In contrast to density-density correlations, vbsa is not very sensitive to the quark 
mass over the small range investigated here and it is roughly a factor of two smaller than 
the experimental value. As discussed above, using a more realistic value for w will reduce 
this deviation. Nevertheless, based on the argument given below we expect the rms radius 
measured from the BSA to be smaller than the experimental value. 

The pion wavefunction on any given spacelike surface can be schematically decomposed 
into its Fock space components made up of quarks and gluons 

k) = ^75^^) + c° |c/75M(ng')) + cl^\dj5uqq{ng)) + ... (5.8) 
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On the lattice the higher Fock states are present because the quark propagator, even in 
the quenched approximation, has an overlap with multi-quark and gluon states due to 
the back and forth propagation across any spacelike surface. The BSA only probes those 
Fock states in which dj^u is traced over spin and color to match the quantum numbers 
of a pion, while the density-density correlation has no such restriction. This distinction is 
crucial because the charge radius gets enlarged if the strong force between the d and u can 
be screened by popping qq pairs from the vacuum. This is shown schematically in Fig. 9 
for the density-density correlation with both the u and d forming a color singlet state 
with a qq pair produced either by a "Z excursion" or by vacuum polarization. Due to the 
restriction of spin and color trace on dj^u only those higher Fock states which factorize, for 
example \dj5uqq(ng)) — )• \dj5u)\qq(ng)) or those that can be rearranged into this form by 
spin and color Fierz, will contribute to the BSA. For this reason the rms radius calculated 
from BSA is expected to be an underestimate of the experimental value. 

Using the same set of propagators and image correction scheme, we have measured 
the pseudoscalar correlation 'R.^^(x) for the pion. The radial distribution x'^'R.^^{x) at 
K = 0.154 is plotted in Fig. 10. The asymptotic behavior of 'R.^^(x) is simpler to analyze 
than 'R.'^'^{x); since the 75 insertions are unchanged under a Wick rotation the exponential 
decay of 'R.^^(x) for large x is controlled by the pion mass. The data in Fig. 10 show that 
a 16^ lattice is too small to observe this asymptotic decay. In Ref. [18] it was shown on 
16"^ lattices at /3 = 5.7, corresponding to a physical volume larger by a factor of two, that 
the asymptotic behavior is consistent with the expected behavior 

The volume integral of p^^ characterizes interactions between quarks. For example in 
the MIT Bag model, Lissia et al. [21] have shown that for hadronic systems constructed 
from noninteracting quarks, i.e. where the composite wavefunction is a direct product 
of S-wave single quark wavefunctions that satisfy the Dirac equation, the volume integral 
J d^ xx'^ p^^ (x) is zero. Deviations of this volume integral from zero result from interactions 
between the quarks. The present data show the requisite node but we cannot estimate the 
integral with sufficient accuracy due to the finite size cutoff, so we cannot test the above 
hypothesis. 
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6. Polarization dependence of the rho meson BSA 

The p meson wavefunction is a linear combination of / = and 1 = 2 orbital angular 
momentum states. On the lattice, where rotational invariance is absent, the BS amplitude 
decomposes into irreducible representations of the cubic group: 

(0Kx)7,A^(0,x)u(0)|K0,j)) 

= ^ (<^..>A,(x) + 8,, (^|| - ^e{x) + (1 - 8,, f^M.x)) ^^'^^ 

where |/9(0,j)) is a state of momentum and polarization j. The functions , <i>E and 
cf)T2 are scalars which multiply tensors transforming under the cubic group as the Ai , 
E and T2 representations (of dimensions 1, 2 and 3 respectively). As is clear from the 
tensor structure, Ai corresponds to the / = state while E and T2 together form the 
decomposition of the 1 = 2 state. 

Lattice calculations therefore allow us to investigate, as a function of the quark mass, 
the relative mixture of / = and 1 = 2 states {Ai versus E or T2), and the restoration of 
rotational symmetry (degeneracy of E and T2) by studying the three cases: (A) i = j and 
X along i (||), (B) i = j and x perpendicular to i (-L), and (C) i ^ j ■ 

At present we only measure the BS amplitude for i = j ^ and the results for cases (A) 
and (B) at «: = 0.154 are shown in Fig. 11. The data show that for large separation the 
fall-off is extremely well fit by an exponential in both cases, with a rate of fall-off given 
by my = 0.226(8) and m± = 0.263(7) respectively. It is interesting to note that the large 
X fall-off is governed by a mass that is roughly mp/2 = 0.23. Such a behavior would be 
expected if the rms radii defined in Eqns. (5.6) and (5.7) were the same. 

From the data shown in Fig. 11 we extract (f>Ai{x) and cf)E{x). The results, shown in 
Fig. 12, are fit to simple hydrogen-like radial wavefunctions 
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(j)Ax{x) = 1.61 e' 

(6.2) 

^e{x) = 0.029 x^ e-'-'''^ . 

These functions give a good fit to (f>Ai{x) for x > 6 and cf>E{x) for x > 2. The results 
at /€ = 0.155 are qualitatively similar. The data for both cf)Ai{x) and cf>E{x) are slightly 
broader, though the difference is smaller than the statistical errors. 
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7. Conclusions 

We show how to define gauge invariant BSA using smeared hnks that have similar 
overlap with meson wavefunctions as those in Coulomb or Landau Gauge. We use the gauge 
invariant BSA to calculate the rms radius for the pion and qualitatively demonstrate the 
Lorentz contraction along the direction of motion for a pion. 

For the rho meson the wavefunction is characterized by three different orbital functions 
(pAi 1 (pE and (f>T2 ■ These correspond to the Ai , E and T2 representations of the cubic group, 
which in the continuum limit become the / = and 1 = 2 representations. We show how 
each of these three functions can be extracted by choosing the polarization and separation 
axis appropriately. A test of restoration of rotational invariance can be made by comparing 
(pE and (pAi ■ With the present data we can only calculate (pAi and (pE and show that these 
are well described by simple / = and 1 = 2 hydrogen-like wavefunctions respectively. 

The rms radius of the pion measured from the BSA is roughly 0.5 — 0.7 of the exper- 
imentally measured value. The large uncertainty is due to the lack of knowledge of the 
position of the center of mass with respect to the valence d and u due to the motion of 
the glue and additional qq pairs. We attribute the remaining discrepancy to the following 
two factors. First, the BSA does not fully probe those parts of the wavefunction in which 
there is a saturation of the color force between the valence quarks and qq pairs in the sea. 
Second, in the quenched approximation, sea quark contributions of the type shown in Fig. 
9b are absent. 

Present data show that the pion's radius extracted from density-density correlations 
has a significant variation with the quark mass. An extrapolation to the physical pion 
mass gives a significantly larger value than the experimental number. This is not a cause 
for concern as rdd is not simply related to the experimentally measured charge radius. 
Furthermore a comparison of data between /3 = 5.7 and 6.0 shows that the deviation 
decreases with the lattice spacing and suggests that taking the continuum limit would give 
a result much closer to the experimental value. Lastly, one needs to evaluate the effects 
of using the quenched approximation, excited state contamination, and the disconnected 
diagrams that have been neglected in this calculation. 
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/3 


n 










6.0 


0.154 


0.365a 
700MeV 


6.8(l)a 
.72(l)fm 


6.5(l)a 
.68(l)fm 


3.01a 
0.32fm 


6.0 


0.155 


0.296a 
560MeV 


7.6(l)a 
.80(l)fm 


7.0(2)a 
.74(2)fm 


3.08a 
0.33fm 


5.7 


0.160 


0.694a 

694MeV 


3.8(l)a 
.76(l)fm 


3.8(l)a 
.76(l)fm 




5.7 


0.164 


0.527a 

527MeV 


4.5(l)a 
.89(l)fm 


4.5(l)a 
.89(l)fm 





Table 1. Results for the pion's rms radius calculated from the BSA amplitude and the 
density-density correlations p^^i^x). The results are given in both lattice units and in 
physical units. We use a = 0.2fm and a = 2/19fm for the lattice scale at /3 = 5.7 and 
6.0 respectively. The estimate for tbsa has 20% uncertainty. The experimental value is 
= 0.636 ± 0.037fm. 
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Fig. la. Gauge invariant BS amplitude for the pion at various time slices from the source. 
The data are obtained after the gauge links have been smeared 6 times as described in the 
text. 
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Fig. lb. Same as Fig. la for the p meson. The data are for the case where the separation 
X is taken along the polarization axis. 
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Fig. 2a. The gauge-invariant BS amplitude at t = 15 as a function of smearing level. 
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Fig. 2b. Same as Fig. 2a but for the p meson. 
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Fig. 3. Exponential fit to the gauge invariant BS amplitude for the pion at large x. The 
rate of fall-off is extracted from a fit to points 6 < x < 12 at each of the two masses. 
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Fig. 4. Comparison of the pion BS amplitude at zero and non-zero momentum. The data 
show Lorentz contraction when x\\p and pa = (0, 0, 27r/16). 
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Fig. 5a. Comparison of Cjr, A^^^ and A^^^ at k = 0.154. The data for at large x have 
not been corrected for contributions from mirror sources. 
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Fig. 5b. Comparison of Ct^-, A^^^ and Ai^\ Rest is same as Fig. 5a. 
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Fig. 6 The three different Wick contractions that contribute to the density-density cor- 
relation. The fiUed circles denote the meson source and sink while the filled rectangles are 
the density insertions. 
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Fig. 7 The two different types of diagrams that give the leading large x behavior for the 
density-density correlation. 
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a) 



b) 




Fig. 9 We show two example of how the color force can be saturated in the the density- 
density correlation. The filled rectangles are the density insertions and the shaded blobs 
are a schematic representation of a color neutral object. The thick line in (a) shows a "Z 
excursion" by the quark propagator. 
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X 

Fig. 10 The fall-off with x of the p^^{x) density-density correlation inside a pion. The 
correlator has a node at x 8.7a. The lattice is not large enough to verify the expected 
asymptotic behavior e~"^^^. 
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Fig. 12. Results for (f>Ai{x) and cf>E{x) using the data shown in Fig. 11. The fits are 
described in Eq. (6.2). 
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